Elastic turbulence in a curvilinear channel flow 
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We report detailed quantitative studies of elastic turbulence in a curvilinear channel flow in a dilute 
polymer solution of high molecular weight polyacrylamide in a high viscosity water-sugar solvent. 
Detailed studies of the average and rms velocity and velocity gradients profiles reveal an emergence 
of the boundary layer associated with the nonuniform distribution of the elastic stresses across the 
channel. The characteristic boundary width is independent of the Weissenberg number Wi and 
proportional to the channel width that follows from our early investigations of the boundary layer 
in elastic turbulence of different flow geometries. The appearance of the characteristic spatial scales 
of the order of the boundary layer width of both velocity and velocity gradient in the correlation 
functions of the velocity and velocity gradient fields in a bulk flow suggests that rare and strong 
parcels of excessive elastic stresses, concentrated in the boundary layer, are ejected into the bulk 
flow similar to jets observed in passive scalar mixing. And finally, the experimental results show 
that one of the main predictions of the theory of elastic turbulence, namely the saturation of the 
normalized rms velocity gradient in the bulk flow of elastic turbulence contradicts to the experimental 
observations both qualitatively and quantitatively in spite of the fact that the theory explains well 
the observed sharp decay of the velocity power spectrum. The experimental findings call for further 
development of theory of elastic turbulence in a bounded container, similar to what was done for a 
passive scalar problem. 

PACS numbers: 47.20.Gv, 47.50.-d, 47.27.-i 



I. INTRODUCTION 

The addition of small amount of long polymer 
molecules into a fluid makes it elastic and capable to 
storingelastic stresses that may strongly alter flow prop- 
erties [l|. First of all, the elastic stresses generated by 
the polymer stretching in the flow becomes the main 
source of nonlinearity in the polymer solution flow at low 
Reynolds numbers Re. As the result, an elastic insta- 
bility shows up, when the elastic energy overcomes the 
dissipation due to polymer relaxation. The ratio of the 
nonlinear elastic term to the linear relaxation is defined 
by the Weissenberg number Wi [l|, which is the main 
control parameter in the problem, and the elastic insta- 
bility occurs in a shear flow with curvilinear trajectories 
at Wi c > 1. Above the purely elastic instability, a path 
to a chaotic flow in a form of irregular flow patterns at 
Wi > Wi c was studied in three flow geometries: Cou- 
ette flow between cylinders, swirling flow between two 
disks, and flow in a curvilinear channel 0, Further 
increase of Wi at vanishingly small Re leads to the most 
remarkable phenomenon discovered recently experimen- 
tally 0| and then studied during the last decade in the 
increasing number of experimental and theoretical 
[i^-fl2| papers, namely "elastic turbulence". It is a spa- 
tially smooth and random in time flow, which is driven 
by strong polymer stretching and resulting elasticity ob- 
served at sufficiently large Wi and at vanishingly small 
Re. Many properties of the elastic turbulence regime, 
namely statistics of velocity and velocity gradient fields, 
spatial and temporal velocity correlation functions and 
power-law behavior of velocity power spectrum, a new 
length scale-the boundary layer width and its properties, 
scaling and statistical properties of torque and pressure 



fluctuations, in a von Karman swirling flow between two 
disks were observed and investigated experimentally and 
only a few of them theoretically and numerically [7] . 

On the other hand, a detailed description of similar 
properties in a channel flow is lacking, though elastic 
turbulence in a curvilinear channel was used to demon- 
strate its effectiveness in mixing in macro- as well as in 
micro-channels 13 . In a 3 mm wide curvilin- 
ear channel the longitudinal and transverse flow velocity 
components were measured by a laser Doppler anemome- 
ter at the bend N = 12 near the middle of the half-ring 
in the middle of the channel at a fixed value of Wi, about 
twice above the onset to the elastic instability. The power 
spectra of the longitudinal and transversal velocity com- 
ponents show a broad region of an algebraic decay /~ 3 ' 3 
in the frequency domain. The rms velocity of fluctuations 
was 0. 09(F) and 0. 04(V) for the longitudinal and trans- 
verse components, respectively, where (V) is the mean 
longitudinal velocity component. Since the power spec- 
tra were measured at a point with more than 10 times 
higher mean than the characteristic fluctuation velocity, 
the Taylor hypothesis can be used to transfer the power- 
law-decay in the frequency domain to the wave number 
k domain as k~ 3 3 , which is very close to the scaling ex- 
ponent found in the von Karman swirling flow between 
two disks in both frequency and wave number domains 
0, 0] ■ Another measurements of the flow velocity were 
carried out in the micro-channel of a similar design as the 
3 mm channel but scaled down 30 times compared with 
the macro-channel version [H, [l3|]. In order to define the 
onset of the elastic instability, the both longitudinal and 
transverse velocity field components were measured by 
microscopic particle image velocimetry, micro-PIV, as a 
function of the pressure drop Ap along the channel. The 



2 



resulting flow resistance, defined as Ap/V, as a function 
of the pressure drop shows sharp but continuous change 
in the dependence that determines the instability onset. 
Similar effect is observed on the plot of Vg ms as well as 
yrms versus [j| [l3[ . Temporal dependence of the lon- 
gitudinal velocity as well as its correlation function were 
also measured together with the velocity correlation time. 
By using 0.2 /im fluorescent particles streamwise vortices 
were visualized by means of horizontal confocal scanning 
microscopy in the middle plane of the channel Il3j . 
In the recent detailed studies of mixing in a curvilinear 
channel of 1 mm 2 cross-section using elastic turbulence, 
some partial characterization of the velocity field neces- 
sary for quantitative characterization of mixing, namely 
both average radial (V r ) and longitudinal (Vg) velocity 
profiles across the channel, dependence of (Vg) on Wi, 
profiles of rms of longitudinal dVg/dr and radial dV r /dr 
velocity gradients, was carried out [l4j . 

The first and currently the only theory of elastic turbu- 
lence was developed right away after the first publication 
of its discovery. The main concern of the theory was to 
explain the key experimental observation in elastic tur- 
bulence, namely the sharp algebraic decay of the velocity 
power spectrum with the scaling exponent 5 between - 
3.3 and -3.6 [2j,|3j. Due to the sharp velocity spectrum 
decay, the velocity and velocity gradient are both deter- 
mined mostly by the integral scale, i.e., the vessel size. 
It means that elastic turbulence is essentially a spatially 
smooth and random in time flow, dominated by strong 
nonlinear interaction of a few large-scale modes. It is 
the same random flow that occurs in hydrodynamic tur- 
bulence below the dissipation scale and called Batchelor 
flow regime [16]. It is resulted from stretching and folding 
of elastic stress field, similar to a passive scalar stretching 
and folding in the Batchelor regime of mixing. A crucial 
difference between elastic turbulence and passive scalar 
mixing is that in the case of elastic turbulence the corre- 
sponding elastic stress field is not passive but reacts back 
on the driving velocity field and in such way stabilizes the 
flow There are two aspects of theory of polymer 

stretching in a flow. First is a description of statistics 
of polymer stretching and a coil-stretch transition in a 
spatially smooth and random in time flow, and second is 
a characterization of the properties of elastic turbulence 
resulting from the polymer stretching. 

The first aspect requires a microscopic approach to the 
problem, which provides a quantitative prediction on the 
coil-stretch transition of a polymer and on a saturation 
of the polymer stretching in a spatially isotropic random 
unbounded flow and a detailed prescription to experi- 
mentally verify it [lU, [l9|. This prediction was tested 
experimentally, and good agreement was found (20l . [2l| . 
The coil-stretch transition has remarkable macroscopic 
consequence on a flow: properties of the polymer solu- 
tion become essentially non- Newtonian and the stretched 
polymers significantly alter the flow due to their back re- 
action. The second aspect, on the other hand, requires a 
macroscopic description of elastic turbulence, which has 



been developed by Lebedev et al [9|, ll7j and is based on 
polymers with linear elasticity and the feedback reaction 
on the flow. The theory of elastic turbulence uses the 
set of equations for the elastic stress tensor and velocity 
fields. Hydrodynamic description of a polymer solution 
flow and of dynamics of elastic stresses for linear poly- 
mers is analogous to that of a small-scale fast dynamo 
in magneto-hydrodynamics (MHD) and also of turbu- 
lent advection of a passive scalar in the Batchelor regime 
[j| Hijj . though some significant differences exist. The 
stretching of the magnetic lines is similar to the polymer 
stretching, and the difference with MHD lies in the re- 
laxation term that replaces the diffusion term in MHD 
description, whereas in the passive scalar advection prob- 
lem the dynamo effect, i.e. feedback reaction on the flow, 
is absent. In all three cases the basic physics is the same, 
rather general and directly related to the classical Batch- 
elor regime of mixing: stretching and folding of the pas- 
sive scalar, magnetic, or stress fields by a random advert- 
ing flow in all three cases. 

Theory of elastic turbulence in an unbounded flow of 
a polymer solution is based on the following assumptions 
[a LLDl ■ (i) A statistically stationary state occurs due to 
the feedback reaction of stretched polymers (or the elas- 
tic stress) on the velocity field that leads to a saturation 
of the elastic stress a p and rms of the velocity gradients 
{dVi/dXj) rms (and so Wii oc = (dVi/dXj) rms X, where A 
is the longest polymer relaxation time) even for polymers 
with linear elasticity [l7| . The saturation value in a 
bulk of elastic turbulence is Wii oc ~ 1 and constant at 
all Wi above the coil-stretch transition. It is the key 
theoretical prediction 0, [TtJ • (ii) Both dissipative terms 
due to viscosity and polymer relaxation, which ap pea r 
in the equation for the dissipation of elastic energy [17| . 
are of the same order, i.e. a p /X ~ ry(Vy) 2 or otherwise 
<J v \jr\ ~ Wif oc , where r\ is the viscosity. Then both as- 
sumptions lead to the following result: the normalized 
elastic stress a p X/r] ~ 1 and also saturates. These two 
important theoretical predictions deserve a stringent ex- 
perimental test. Indeed, the value of elastic stresses was 
measured in the recent experiment in microscopic ar- 
rangement of a swirling flow [l5[ and strong discrepancy 
was identified. The first prediction on the saturation of 
Wii oc in a bulk flow of elastic turbulence was examined in 
a swirling macroscopic flow [?J and quantitative disagree- 
ment was found. On the other hand, the further theoret- 
ical analysis leads to a power-like decaying spectrum for 
the elastic stresses and for the velocity field fluctuations 
with the exponent \S | > 3 in a good accord with the ex- 
perimental results @, [1] ■ The close value of the exponent 
in the velocity power spectra decay was also obtained in 
the numerical simulations of elastic turbulence based on 
the Kolmogorov shear flow of a dilute polymer solution 
described by the Oldroyd-B model [IS Ell- 
in this paper we provide a complete characterization of 
the channel flow in the elastic turbulence regime, present 
statistics, correlations and scaling of velocity and veloc- 
ity gradient fields in the bulk as well as in the vicinity 



of the wall, and also properties of the velocity gradie 
boundary layer that provide us the possibility also to fu 
ther test the theoretical prediction about the saturatk 
of Wiioc. 

The paper is organized as follows. In Sec. II, the e 
perimental setup is described in details. In Sec. Ill tl 
experimental results are presented. Here we first descril 
in Sec. IIIA flow structure and experimental determin 
tion of the elastic instability threshold. Then in Se 
IIIB various velocity and velocity gradient profiles ai 
determination of the boundary layer widths of veloci 
and velocity gradient and its functional dependence a 
described. Temporal and spatial correlation functions 
velocity and velocity gradients and correlation times ai 
lengths as a function of Wi are presented in Sec. III< 
And finally, statistics of velocity and velocity gradien 
and structure function scalings are discussed in Sec. III! 
Discussion of the experimental results is given in Sec. I 
and conclusions are presented in Sec. V. 



II. EXPERIMENTAL SET-UP AND 
PROCEDURE 

The experiments were conducted in a curvilinear cha 
nel of about 1 mm 2 square cross-section machined out 
plexiglass (lucite). The width size was chosen to redu 
by an order of magnitude the amount of a working flu 
(polymer solution) used in the first experiment with 
channel width of 3 mm [j| H[ and to increase spatial re 
olution of a velocity field in peripheral regions compared 
with the second experiment, where a micro-channel of 
100 /j,to 2 square cross-section was used [!, HH. The chan- 
nel used in the current experiment contained 63 pairs of 
identical smoothly connected half-rings (bends or units) 
with inner and outer radii of Ri — 1 mm and Rq = 2 mm, 
respectively, and sufficiently high gap ratio d/Rt = 1, 
the same as in the previous experiments [Tij , which was 
intended to facilitate the onset of the elastic instability 
at sufficiently low Wi. Here d = Rq — Ri is the chan- 
nel width. The precise dimensions of the channel were 
0.95 mm width at the midplane of the channel, where 
measurements of mixing were performed (with ±0.05/im 
differences in the width at the top and bottom of the 
channel), and 1.025 mm depth. Thus the entire chan- 
nel length was approximately 59.4 cm measured along a 
channel midplane (see Fig. 1). The channel main body 
and the lid sealed with O-ring were squeezed between two 
stainless steel plates to seal the channel against leaks and 
to preserve a flatness along the channel. 

The pulse Nd-YAG laser of the 532 nm wavelength 
(New Wave Research Ltd.), which produces pulses of 
power 30 mJ/5ns with time delay between the pulses of 
10 /us, was used to conduct particle image velocimetry 
(PIV) measurements. Time differences between two con- 
secutive images was chosen depending on Wi between 
0.51 and 1.84 msec, and it was synchronized between 
laser pulses and camera via the control units of the PIV 
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Figure 1: The schematic drawing of the experimental setup 



laser system. For visualization of velocity field measure- 
ments, red fluorescent particles (Duke Scientific Ltd.) of 
2 /xm at concentration of 150 ppm were used. As in the 
previous setups 0, |3, [HJ, the channel was illuminated 
from a side by laser beam transformed by an appropriate 
optical setup containing two cylindrical lenses to a broad 
light sheet of a thickness about ~ 50 /im in the obser- 
vation region (see Fig. 1). Thinner light sheet causes a 
permanent damage to plexiglass. The laser beam sheet 
produced a thin cut in the 3D flow, parallel to the top and 
bottom of the channel at its midplane. Fluorescent light 
emitted by the particles in the direction perpendicular 
to the beam plane was detected by a charge coupled de- 
vice (CCD) via a microscope objective EF 4.0/0.12 (Leitz 
Wetzlar) mounted on the plastic tube (guide) . The CCD 
camera Grasshopper, model Gras-1455M-C (Grey Point 
Research) of 16 bits grey scale resolution and with spa- 
tial resolution of 1280 x 768 pixels at up to 30 frames/sec 
rate was used. Together with the microscope objective, it 
provided a PIV spatial resolution of 0.576 /im/pixel. The 
size of the PIV image was 1280 x 384 pixels. The window 
size to get velocity vectors was taken 32 x 8 pixels that 
corresponded to 18.4 x 4.6 nm 2 . Total 130 velocity vec- 
tors in the transverse direction to the flow were obtained 
from PIV. Thus the velocity gradients were calculated on 
dr = 4.6/im. The images were taken at the rate 10 fps. 
Since each two consecutive images were transformed into 
one velocity field, the final rate was 5Hz of the velocity 
field sampling. Total 1000 velocity field were used for 
averaging. 

The velocity profile in the channel flow of the Newto- 
nian solvent measured by PIV has been used for calibra- 
tion. Figure 2 presents the results of longitudinal veloc- 
ity profiles Vg(r/d) across the channel in a laminar flow 
of the solvent for various pressure drops. (We use in the 
channel middle plane as the coordinates z and r, as longi- 
tudinal and transversal to the channel walls, respectively, 
and Vg and V r as the longitudinal and radial components 
of the velocity field, respectively.) We visualize just a 
part of the channel cross-section to increase a spatial res- 
olution, since we later on use only the velocity maximum 
values and the velocities in a peripheral region close to 
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Figure 2: Longitudinal velocity profiles across a channel of a 
solvent laminar flow in a curvilinear channel at various pres- 
sure drops (starting from small values at bottom to top; r — 
corresponds to the inner channel wall). 




Figure 3: Maximum of the velocity profile versus pressure 
drop along the channel AP for the solvent of water with 65% 
of sucrose in a laminar channel flow. 



the wall. Due to curvature in serpentine geometry of the 
channel, the profiles of the longitudinal velocity are not 
symmetric relatively to the middle plane in contrast to a 
straight channel. Figure 3 shows a linear dependence of 
the maximum values of the velocity profiles V g max (or the 
discharge) in a laminar flow of the solvent as a function 
the pressure drop AP along the channel, which is used 
for calibration. 

As a working fluid, a 65% sucrose-water solvent with 
addition of 1% NaCl, 250 ppm of NaN s , and 80 ppm by 
weight high molecular weight polyacrylamide (PAAm) 
(M w = 18 Mda, supplied by Polysciences), was used. 
It was prepared from a master water solution contained 
3000 ppm of PAAm, 1% NaCl, 250 ppm of NaN 3 , and 
3% of iso-propanol [jjj . The viscosity of the solvent 



r/ s = 113.8 mPa ■ s and the polymer solution r\ = 137.6 
mPa ■ s were measured at 22 ± 0.5 °C, the same temper- 
ature kept during the experiment. The longest polymer 
relaxation time measured by the stress relaxation method 
was A = 11.5 sec and was found to be independent of the 
shear rate [l5[ (it differs from values used in Ref . [l3[ , 
where the polymer relaxation time was measured by the 
small amplitude oscillation method). The solution den- 
sity was p = 1.303 ± 0.03 g cm -3 . The inflow of a poly- 
mer solution into the channel was generated via two in- 
lets by a compressed nitrogen gas at pressures between 
0.35 and 8.5 psi depending on Wi. The gas pressure was 
regulated and measured via a regulated pressure gauge 
with a precision of ±5 Pa. The compressed gas was fed 
into two plastic cylinders containing a working polymer 
solution. The cylinders were connected to the channel 
inlets by Tygon tubes with inner diameter of 1.0 mm, 
sufficiently large to prevent a possibility of an elastic in- 
stability already in the tubes. The Weissenberg numbers 
Wi = (2Vg nax / d)X reached in the experiment by a such 
arrangement were in the range from to 950. The max- 
imum Re = V e max dp / r\ was ~ 0.54, so that the inertial 
effects in the flow were negligible, as was already shown 
in the earlier experiments [3j-|5j, Il3| . 



III. PROPERTIES OF ELASTIC TURBULENCE 

A. Flow structure and experimental determination 
of the elastic instability threshold 

We searched experimentally for the onset of the elastic 
instability in our channel with the aspect ratio dj ' Ri = 1. 
In order to figure out the range of the elastic turbu- 
lence region, we plot the values of (Vg) at r/d = 0.5 
and V g rms as a function of Wi in Figs. 4 and 5. From 
the both plots one concludes that the elastic instability 
onset corresponds to the change in the slope that occurs 
at Wi c ~ 200 on the both plots. Similar change in the 
slope takes place in the dependence of V" ns on Wi at 
about the same value of Wi c though the data is noisier 
(Fig. 6). An additional change in functional behavior of 
V e rms and V r rms versus Wi occurs at Wi w 350 - 400, 
where the transition region to elastic turbulence is ended. 
In all three plots in Figures 4, 5, and 6 the arrows indi- 
cate the onset of the elastic instability and the end of the 
transitional region towards elastic turbulence. Thus, the 
range of developed elastic turbulence can be considered 
for Wi from about 350 up to 950. 

Using the results of the PIV measurements one cal- 
culates the various components of velocity gradient field 
as a function of r at different values of Wi. Figure 7 
presents the data for Wii oc = (dVe/dr) rms X versus Wi 
averaged over r/d from 0.2 till 0.5 in the bulk flow. In the 
range of elastic turbulence for Wi from 350 till 950 the 
data in log-linear presentation in Fig. 7 are fitted by a 
linear fit, so Wi loc = 82.2 exp(Wi/ 919. 2). Figure 8 show 
(dV r /dr) rms averaged over r/d from 0.2 till 0.5 in the 
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Figure 4: The mean longitudinal velocity at r/d = 0.5 as a 
function of Wi. Arrows indicate the elastic instability onset 
and the ending of a transitional regime to elastic turbulence, 
respectively. 



Figure 6: Root mean square of the transversal velocity V" ns 
versus Wi. Arrows indicate the elastic instability onset and 
the ending of a transitional regime to elastic turbulence, re- 
spectively. 
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Figure 5: Root mean square of the longitudinal velocity Vg~ ma 
versus Wi. Arrows indicate the elastic instability onset and 
the ending of a transitional regime to elastic turbulence, re- 
spectively. 




1000 



Figure 7: The local Weissenberg number Wii oc = 
(dVe/dr) rm3 \ averaged over r/d from 0.2 till 0.5 as a function 
of Wi. Arrows indicate the elastic instability onset and the 
ending of a transitional regime to elastic turbulence, respec- 
tively. 



bulk flow region as a function of Wi. Here the data can 
be also fitted linearly by (dV r /dr) rms = 5.64 + 0M5Wi 
in the elastic turbulence range. All PIV data were aver- 
aged over 1000 velocity fields to get both average velocity 
and its rms values as well as the velocity gradient com- 
ponents. 



B. Velocity and velocity gradient profiles and the 
boundary layer problem 

As it is demonstrated above, a channel flow of a poly- 
mer solution at sufficiently large Wi above the elastic 



instability threshold Wi c is chaotic with large velocity 
fluctuations. Figure 9 shows the profiles of the average 
longitudinal velocity (Vo)(r/d) across the channel for var- 
ious values of Wi from below the elastic instability and 
up to the highest values. A boundary layer character- 
ized by a sharp drop of (Vg)(r/d) at about r/d < 0.05 is 
clearly seen at Wi above 192. A full profile of (V$)(r/d) 
across the channel at Wi — 951 is presented in Fig. 10, 
where both boundary layers are clearly identified. 

From the same velocity measurements the profiles of 
the average transversal velocity (V r ), and the rms fluc- 
tuations of the longitudinal Vg ms and transversal V™ s 
components of the velocity for various Wi are also ob- 
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Figure 8: (dV r /dr) rms averaged over r/d from 0.2 till 0.5 
versus Wi. Solid line is a linear fit to the data (dV r / dr) rrns = 
(5.64 ± 0.15) + (0.055 ± 0.002)Wi in the elastic turbulence 
regime from Wi = 350 till 950. 
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Figure 9: Mean longitudinal velocity profiles across the chan- 
nel of a polymer solution flow for various Wi (starting from 
small values at bottom to top; r = corresponds to the inner 
channel wall). 



tained (see Figs. 11,12,13). 

First, by comparing Figs. 9 and 11 one finds that 
(V r ) is close to zero in the limit of error bars across the 
channel and (Vg) exceeds it by more than two orders of 
magnitude. Whereas rms fluctuations of the longitudi- 
nal and transversal components are of the same order 
of magnitude that follows from Figs. 12 and 13. Fig- 
ure 14 shows (dVff/dr) rms versus r/d at various Wi in 
the whole range of its variation. Analogously, the com- 
ponent (dV r / dr) rms as a function of r/d is plotted for 
various Wi in Fig. 15. As before, all PIV data were 
averaged over 1000 velocity fields to get both average ve- 
locity and its rms values as well as the velocity gradient 



Figure 10: The entire mean longitudinal velocity profile across 
the channel of a polymer solution flow for Wi = 951. 
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Figure 11: Mean transversal velocity profiles across the chan- 
nel of a polymer solution flow for various Wi. 



components. 

Second, we present in Figs. 16, 17 time dependencies 
of Vg(t) and V r (t) for fixed value of r/d — 0.5 in the 
bulk flow at several values of Wi and on time interval 
much larger than A. In the elastic turbulence regime at 
Wi > 350 strong fluctuations in the both components of 
the velocity are observed. 

Clearly identified boundary layer regions near the wall 
in the both average velocity and rms of fluctuations of 
both components of the velocity and velocity gradient 
profiles across the channel in Figs. 9,10,14,15 can be 
better study after subtraction from each velocity profile 
a linear part of the profile found in the bulk region with 
the slope d{Vg)/dr (see Figs. 9 and 10) and then scaling 
its value to unity. As the result, all longitudinal velocity 
profiles at Wi in the range of elastic turbulence above 350 
collapse on a single curve demonstrating independence 
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Figure 12: RMS fluctuations of the longitudinal velocity pro- 
files across the channel of a polymer solution flow for various 
Wi (starting from small values at bottom to top; r — cor- 
responds to the inner channel wall). 



Figure 14: (dVe / dr) Tma versus r/d for various Wi from bot- 
tom to top. 
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Figure 13: Rms fluctuations of transversal velocity profiles 
across the channel of a polymer solution flow for various Wi 
(starting from small values at bottom to top; r — corre- 
sponds to the inner channel wall). 
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Figure 15: (dV r / dr) rma versus r/d for various Wi from bot- 
tom to top. 



C. Temporal and spatial correlation functions of 
both velocity components and their gradients. 



of the boundary layer width on Wi (Figs. 18,19). This 
property of the boundary layer agrees with that observed 
in the swirling flow Q • The boundary layer width is also 
independent of Wi\ oc (Fig. 19, inset). 

And finally, we present in Fig. 20 the dependencies 
of (dVg / dr) rms in the bulk and and its peak values in 
the boundary layer regions as a function of Wi taken 
from the plots in Figs. 7 and 14. The peak values of 
(dV$/dr) rms in the boundary layer exceed those in the 
bulk up two times only compared to at least an order of 
magnitude in the swirling flow 

0, HI, probably due to 
the limiting spatial resolution of micro-PIV. The latter 
smears out the peak value of (dVg / dr) rms . 



We studied both temporal and spatial (across the 
channel) correlation functions of the both velocity com- 
ponents. Figure 21 shows temporal correlation func- 
tions of the longitudinal and transversal velocity compo- 
nents C{t) = (5V{t)SV{0))/({6V{0)) 2 ), where 6V e {r) = 
V B (r) - (Vg), and SV r (r) = V r (r) - (V r ) at several Wi 
values in the elastic turbulence regime taken at the bend 
N = 42. The corresponding correlation time calculated 
as r corr = J tC T (t)dt/ J C T (t)dt for the both velocity 
components as a function of Wi are presented in Fig. 
22. As can be seen in Fig. 22, the correlation time for 
8Vg(T) is considerably larger than for for SV r {T). 

The spatial correlation functions C(r/d) = 
(SV{r/d)SV(0))/({6V(0)) 2 } and the corresponding 
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Figure 16: Time series of the longitudinal velocity at r/d 
0.5 for various Wi. 



Figure 18: Normalized and scaled down (Ve)(r/d) velocity 
profiles at various Wi. 
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Figure 17: Time series of the transversal velocity at r/d — 0.5 
for various Wi. 
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Figure 19: Velocity boundary layer width lu versus Wi. Inset: 
Velocity boundary layer width Iti versus Wii oc - 



correlation lengths for the both velocity components cal- 
culated as l corr /d — J (r/d)C(r/d)dr/d/ JC(r/d)dr/d 
are shown for several Wi in Figs. 23 and 24. Here 
5V e (r/d) = Ve(r/d)-(V e ) and SV r (r/d) = V r (r/d)-{V r ). 
The correlation lengths for both components are rather 
close and independent of Wi. Using the Taylor hypoth- 
esis one can compare the correlation lengths along the 
flow with those across the flow. As follows from the 
data in Fig. 22, the correlation time for the longitudinal 
component of velocity is about T corr /X ~ 4 that gives 
for, let say, Wi — 600 the correlation length along the 
flow T/d = T corr (Vg) ~ 370 that is about 2000 times 
larger than l cor r/d. Similarly, the correlation functions 
Cg(r/d) and C r (r/d) as well as the corresponding corre- 
lation lengths le,corr/d and l r . CO rr/d of the both velocity 
gradients dVg/dr and dV r /dr are also calculated, and 
the results for different Wi were presented in Figs. 25 



and 26. First, the correlation functions have minimum 
at r/d w 0.06, and second, the correlation lengths are up 
to an order of magnitude shorter than for the velocities 
(see Fig. 26 versus Fig. 24) and grow with Wi. 



D. Statistics of velocity gradients and their 
structure functions. 

Further verification of the scaling laws in elastic tur- 
bulence in the channel flow comes from the statistical 
analysis of the velocity field. We point out that the sta- 
tistical analysis of the longitudinal velocity component 
in spatial and in particular temporal domains show large 
scatter due to insufficient data that makes the analysis 
unreliable. On the other hand, the same analysis in a 
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Figure 20: The peak values of (dV e /dr) rms in the bulk and 
boundary layer versus Wi. 



Figure 22: Normalized correlation times Tg^ corr /\ and 
T r ,corr/^ versus Wi for both Ve and V r , respectively. 




r/d 



Figure 21: Temporal correlation functions Cb(j) of V$(t) and 
temporal correlation functions CV(t) of V r (t), respectively, at 
several Wi above Wi c taken at the bend N = 42. 



Figure 23: Spatial correlation functions Cg(r/d) of Ve(r/d) 
and spatial correlation functions C r (r/d) of V r {r/d), respec- 
tively, at several Wi above Wi c taken at the bend TV = 42. 



spatial domain of both components of the velocity gradi- 
ents exhibits much better results. 

First, we conducted the statistical analysis of the spa- 
tial increments of the radial gradients of the longitu- 
dinal velocity 6dVeJdr{5r/d) = [dV e /dr(r/d + Sr/d) - 
dVg / d{r / d)]/ SdVg / dr{5r / d) rms in a wide range of the 
spatial scales from 9.2 to 46 /im with the step 4.6 /im 
at Wi = 951 for the bend N — 42. The correspond- 
ing PDFs of the spatial increments of the normalized 

longitudinal velocity gradients 5(dV$ / 'dr)(Sr / d) have a 
small Gaussian cap and show well-pronounced exponen- 
tial tails, clear scale invariance and symmetry with small 
scatter in spite of low statistics (see Fig. 27). Further 
analysis can be done in equivalent way either by direct 
calculations of the structure functions or by calculations 



of the moments of PDFs. The corresponding second mo- 
ments of PDFs S2,e(Sr/d) on 8r/d in log- log coordinates 
for several values of Wi show the scaling region in Sr/d 
between 0.01 and 0.05 in Fig. 28. The structure func- 
tions of the higher even orders up to p = 8 S p ,g(Sr/d) as 
a function of Sr/d are plotted in log-log coordinates at 
Wi = 951 in Fig. 29. Due to symmetrical shape of PDFs 
odd moments are zero. The power dependence of the 
structure functions (or moments) S Pl $(6r/d) ~ (6r 
is found in the range of Sr/d between 0.004 and 0.05. 
The plot in Fig. 30 demonstrates independence of the 
scaling exponents £2,0 and of Wi in the whole range 
of elastic turbulence. The dependence of £ Pi g which are 
surprisingly close to the linear scaling with £ p g = 0.75p 
(see Fig. 31). The latter is very different from passive 
scalar behavior [14]. Analogous analalysis was conducted 



10 




400 500 600 700 800 900 1000 



Wi 




1000 



Figure 24: Normalized correlation length le,corr/d and 
versus Wi for both V$ and V r , respectively. 



Figure 26: Normalized correlation length le. CO rr/d and 
lr,corr/d versus Wi for both dVe/dr and dV r /dr, respectively. 
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Figure 25: Spatial correlation functions Ce(r/d) of dVe/dr 
and spatial correlation functions C r (r/d) of dV r /dr, respec- 
tively, at several Wi above Wi c taken at the bend N = 42. 
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Figure 27: PDFs of the spatial increments of the normal- 
ized longitudinal velocity radial gradient 8(dVg/dr(8r /<£)) = 
[dV 9 /dr{r/d + 5r/d) - dV e /dr(r/d)]/S{dV e /dr(Sr/d)) rms at 
different length scales (from 9.2 to 46 /xm with the step 4.6 
fj,m) at Wi = 951 based on its spatial velocity field measure- 
ments for the bend TV = 42. 



also for the radial gradients of the transversal velocity 

S(dV r /dr)(Sr/d) and the obtained results were very sim- 
ilar to the velocity gradient of the longitudinal compo- 
nent. The corresponding PDFs also exhibit similar fea- 
tures: small Gaussian cap, exponential tails though not 
so clean, scale invariance and symmetry (see Fig. 32). As 
the result, the second S2,r(Sr/d) and higher order even 
S Pi r(Sr/d) moments (or structure functions) up to p = 8 
show scaling in the same range of scales (Figs. 33, 34, 
35) with the scaling exponents £ p , r and Wi (see Fig. 35) 
and p (see Fig. 36) dependence close to those found for 

6(dVeJdr)(Sr/d). 



IV. DISCUSSION 

Let us summarized first the main observations reported 
above. 

(i) Well-defined threshold of the elastic instability in 
a curvilinear channel flow is identified by the average 
longitudinal velocity and rms longitudinal and transver- 
sal velocity fluctuations (Figs. 4-6). The transition to 
elastic turbulence regime is determined from Vg ms and 
yrms ag we jj ag (^gVe/dr) rms , or Wii oc , dependencies on 
Wi (Figs. 5-7). In the latter case, an exponential de- 
pendence of Wii oc on Wi above the transition to elas- 
tic turbulence is observed (Fig. 7). (ii) The profiles of 
the average longitudinal velocity are altered drastically in 
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Figure 28: Second moments S2,e(Sr/d) of PDFs of the longi- 
tudinal velocity gradient increments 8(dVg/dr(Sr/d)) versus 
Sr/d for several values of Wi (in log- log coordinates). 



Figure 30: Scaling exponents of the second S2,e(Sr/d) and 
fourth S4.e(Sr/d) moments of the longitudinal velocity gradi- 
ent increments S(dVe / dr(Sr / d)) for various Wi. 
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Figure 29: Structure functions S P] e(5r/d) of the longitudinal 

velocity gradient spatial increments S(dVg/dr(8r/d)) up to 
p — 8 (only even) for Wi = 951 (in log- log coordinates). 



Figure 31: Scaling exponents £ Pt e of the structure func- 
tions of the longitudinal velocity gradients spatial increments 

8(dVe/dr(Sr/d)) versus p for various Wi. 



elastic turbulence compared to a laminar flow (Fig. 9). 
In the elastic turbulence regime the characteristic bound- 
ary layer is clearly identified near the wall (Fig. 9 and 
10). Being normalized and rescaled, all velocity profiles 
are collapsed on one curve with a horizontal flat part in 
the bulk and sharp change near the wall, so that all pro- 
files have the same boundary layer width independent of 
Wi (Figs. 18 and 19). (hi) On the other hand, a profile 
of (dVg /dr) rms (Fig. 14) shows a peak much closer to the 
wall, inside the velocity boundary layer discussed above, 
which location near the channel wall is also independent 
of Wi. As we discussed in our early papers on elastic tur- 
bulence in a swirling flow 0,13, the peak in (dV e /dr) rms 
means also the maximum in elastic stresses, and so the 
boundary layer is defined by the non-uniform spatial dis- 
tribution of the elastic stresses across the channel, (iv) 



The correlation times determined from the temporal cor- 
relation functions for both velocity components differ up 
to 5 times. The correlation time for the transversal ve- 
locity component 7> iCorr is of the order of the polymer 
relaxation time, whereas the correlation time of the longi- 
tudinal velocity component Tg tCorr is several times larger 
(Figs. 21 and 22). (v) The correlation lengths found 
from the spatial correlation functions of both velocity 
components are of about l CO rr/d ~ 0.18 for both veloc- 
ity components in the elastic turbulence regime (Figs. 
23 and 24). On the other hand, the correlation lengths 
obtained from the spatial correlation functions of the ra- 
dial gradients of the longitudinal and transversal velocity 
components are about 6 times smaller than for the veloc- 
ity components in the whole range of elastic turbulence 
(Figs. 25 and 26). (vi) PDFs of the spatial increments 
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Figure 32: PDFs of the spatial increments of the normal- 
ized transversal velocity radial gradient S(dV r /dr(8r /d)) = 
[dV r /dr{r/d + 8r/d) - dV r /dr(r /d)]/8(dV r /dr(8r/d)) rms at 
different length scales from 9.2 to 46 /im with the step 4.6 fim 
at Wi = 951 based on its spatial velocity field measurements 
for the bend N = 42. 
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Figure 34: Structure functions S p , r (r/d) of the transversal 

velocity gradient increments S(dV r /dr(8r/d)) up to p = 6 
(only even) for Wi = 951 (in log- log coordinates). 
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Figure 33: Second moments S2,r(Sr/d) of PDFs of the 

transversal velocity gradient increments 5(dV r / dr(5r / d)) ver- 
sus 8r/d for several values of Wi (in log-log coordinates). 



of the radial gradients of the longitudinal and transver- 
sal velocity components in a wide range of the length 
scales up to the boundary layer width demonstrate the 
scale invariance and exponential tails. It reminds very 
much the properties of the PDFs of passive scalar in the 
Batchelor regime [14j. On the other hand, the second and 
higher order even moments of the PDFs, contrary to the 
passive scalar logarithmic dependence on spatial scales, 
show an algebraic increase with spatial increment with 
the scaling exponents Cp,e & n d Cp,r, which dependence on 
p mildly deviates from linear dependence. 

There are several important messages, which follow 
from the observations summarized above, (i) The elastic 
instability transition in a curvilinear channel is continu- 




1000 



Figure 35: Scaling exponents of the second S2, r (8r/d) and 
fourth S4 t r(8r/d) moments of the transversal velocity gradient 

increments 8(dV r /dr(Sr/d)) for various Wi. 



ous one (forward bifurcation) as already found in early 
experiments [3, [H, El, m contrast to those observed in 
Couette- Taylor flow [4| and swirling flow between two 
disks [1, 0|- (ii) Contrary to the predictions, Wii oc grows 
with Wi in the elastic turbulence regime and its value 
exceeds the theoretically predicted an order of unity by 
more than two orders of magnitude, (iii) The profiles of 
the average longitudinal component of the velocity reveal 
the boundary layer, which width is independent of Wi in 
the whole range of elastic turbulence and being scaled by 
the channel width also independent on the characteristic 
size of the system. Both existence of the boundary layer 
and independence of its width on the control parameters 
of the flow are rather surprising features, in particular 
taking into account that the largest scale of the flow is 
smaller than the dissipation scale. This fact was already 
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Figure 36: Scaling exponents £ PtT of the structure functions of 

the transversal velocity gradients increments 5(dV r / dr(Sr / d)) 
versus p for various Wi. 



twice larger than the velocity boundary layer lu- Since 
elastic turbulence is a smooth random flow, where only 
a few large scale modes dominate the dynamics, one ex- 
pects that l corr should be of the order of d. On the other 
hand, the correlation length defined from the correlation 
function of the velocity gradients l CO rr is about an order 
of magnitude smaller than l cor r- This fact points out on 
the relation of l cor r with the characteristic spatial scale 
corresponding to the peak location of (dVg/dr) rms near 
the wall inside the velocity boundary layer width. Then 
due to eruption of jets of elastic stresses this character- 
istic scale is observed in the bulk, (vi) The same range 
of spatial scales are found in the scaling region of the 
structure functions of the velocity gradients S p that once 
more indicates a possible influence of jets protruding into 
the bulk of the flow. 



V. CONCLUSIONS 



reported in Refs. [y, [7] for the swirling flow, (iv) The 
boundary layer observed in the average velocity profiles is 
a reflection of a nonuniform distribution of the rms of the 
radial gradient of the longitudinal velocity. Indeed, the 
latter profiles exhibit more intricate behavior with sharp 
peaks inside the velocity boundary layer, which locations 
are independent of Wi but the peak values grow with 
Wi. Since (dVe/dr) rms controls the degree of polymer 
stretching in a random flow and in this way the elastic 
stress, one concludes that similar nonuniform distribu- 
tion of the elastic stress can be expected near the wall in 
elastic turbulence. It is a subject for future experiments. 
The expected nonuniform distribution of the elastic stress 
reminds a nonuniform distribution of passive scalar in a 
bounded channel flow |13l . Il4j . l23l | with rare and strong 
ejection of jets occurring in the diffusion boundary layer 
and protruding into the peripheral re gion and even fur- 
ther into the bulk of the channel flow [H [H, EH- These 
jets are considerably alter mixing significantly reducing 
its efficiency. One can expect an emergence of similar 
jets of the elastic stresses, concentrated near the wall in 
the boundary layer, and ejected into the bulk. In such a 
way the elastic stresses are introduced into the flow, (v) 
Another characteristic spatial scale in the flow is the cor- 
relation length of the velocity field l corr , which is about 



To conclude, the experimental results show that one of 
the main predictions of the theory of elastic turbulence, 
namely the saturation of Wii oc in the bulk flow of elastic 
turbulence contradicts to the experimental observations 
both qualitatively and quantitatively in spite of the fact 
that the theory explains well the observed sharp decay of 
the velocity power spectrum 0, [l7| . The nonuniform dis- 
tribution of (dVg I dr) rms across the channel points out on 
the nonuniform distribution of elastic stresses. The lat- 
ter may lead to the rare and strong eruption of the jets 
of elastic stresses from boundary layer into the bulk and 
in this way to introduce small spatial scales into the bulk 
flow. The existence of the resulting velocity boundary 
layer width and its peculiar properties require a proper 
theoretical description. Thus the experimental findings 
call for further development of theory of elastic turbu- 
lence in a bounded container, similar to what was done 



for a passive scalar problem [23 1. 



Acknowledgments 

This work is supported by grants from Israel Science 
Foundation and Lower Saxony Ministry of Science and 
Culture Cooperation Grant. 



[1] R. B. Bird et al. Dynamics of Polymer Liquids, (Wiley, 

New York, 1987). 
[2] A. Groisman and V. Steinberg, Nature 405, 53 (2000). 
[3] A. Groisman and V. Steinberg, Nature 410, 905 (2001). 
[4] A. Groisman and V. Steinberg,New J. Phys. 6, 29 (2004). 
[5] T. Burghelea, E. Segre, I. Bar- Joseph, A. Groisman, V. 

Steinberg,Phys. Rev. E69 066305 (2004). 
[6] T. Burghelea, E. Segre, V. Steinberg, Phys. Rev. Lett. 

96, 214502 (2006). 



[7] T. Burghelea, E. Segre, V. Steinberg,Phys. Fluids 19, 
053104 (2007). 

[8] Y. Jun and V. Steinberg, Phys. Rev. Lett. 102, 124503 
(2009). 

[9] A. Fouxon, V. Lebedev, Phys. Fluids 15, 2060 (2003). 
[10] S. Berti, A. Bistagnino, G. Boffetta, A. Celani, S. Musa- 

chio, Phys. Rev. E 77, 055306(r) (2008). 
[11] S. Berti and G. Boffetta, Phys. Rev. E 82, 036314 (2010). 
[12] H. Yatou, Phys. Rev. E 82, 036310 (2010). 



14 



[13] T. Burghelea, E. Segre, and V. Steinberg, Phys. Rev. 

Lett. 92, 164501 (2004). 
[14] Y. Jun and V. Steinberg, Phys. Fluids 22, 123101 (2010). 
[15] Y. Liu, Y. Jun, V. Steinberg, J. Rheology 53, 1069 

(2009). 

[16] G. K. Batchelor, J. Fluid Mech 5, 113 (1959). 
[17] E. Balkovsky, A. Fouxon, V. Lebedev, Phys. Rev. E 64, 
056301 (2001). 

[18] E. Balkovsky, A. Fouxon, V. Lebedev, Phys. Rev. Lett. 

84, 4765 (2000). 
[19] M. Chertkov, Phys. Rev. Lett. 84, 4761 (2000). 



[20] S. Gerashchenko, C. Chevallard, V. Steinberg, Europhys. 

Lett. 71, 221 (2005). 
[21] Y. Liu and V. Steinberg, Europhys. Lett. 90, 44002 

(2010). 

[22] Y. Jun and V. Steinberg, to be published 
[23] M. Chertkov and V. Lebedev, Phys. Rev. Lett. 90, 
034501 (2003). 

[24] A. Chernykh and V. Lebedev, JETP Lett. 87, 682 (2008). 
[25] H. Salman, P. H. Haynes, Phys. Fluids 19, 067101 (2007). 



